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OPTIMIZATION OF THE FIN GEOJGTHY AND INTIE-FIN 
DISTANCE FCH FODCEB CONVECTION LAMHOi FLO'7 

In this dissertation the optimum fin geometry and 
inter-fin distance have been obtained using recently developed 
optimization techniques. The volume of fins is sought to be 
mininun subjected to the following constraints : 

1. Increased heat transfer in comparison to that cf 
bare tube. 

2. The boundary layer interference. 

3. Biot number for optinun utility of the fin surfaces^ 

4. Number of fins. 

5. Geometrical constraints. 

Available algorithms cf the constrained minimization have 
been applied to the heat transfer problem. 

The electrical network analogy is used to find the tempe- 
rature drop in the fin. The heat transfer rate has been calculated 
by using the electrical analogy method. 

For the optimum fin, with minimum volume, heat transfer 
rate is increased to 40-50 tines that of the hare tube. 

The temperature drop and the heat transfer rate computed 
by the network analogy method are found to be in excellent agree- 
ment with those of numerically computed values. 



CHAPTER 1 


iitf.aducti:." 

1,1 IHTE52UGTICN 

Any heat transfer surface which has metallic projections 
meant to increase the surface area available for heat transfer, is 
called a finned surface . The metallic projections are referred to 
as Ti ns « tube provided with fins is called a finned tube . 

Different geometries, such as spines of triangular, 
circular, rectangular (with constant or variable thickness), 
circular wedge shapes etc. can be used for fins. Annular and 
constant thickness rectangular fins are simple to manufacture. 

Fjns may be transverse or longitudinal. 

Fins may be welded to the primary heat transfering surface 
or integrally cast. Sometimes metal ribbons are peened to tubes. 
Circular fins may be shrun k fitted along the tubes. 

An improved grooved and peened fin tube known as key fin. 
has been developed by the Hunt Heat Exchanger Company jl^. It gives 
a positive bond between the extended and priaaiy surfaces. Soft 
brazing is most common to ensure contact. 

It is argued that the substitution of porous metal for 
fins quadruples the efficiency of heat exchangers and evaporators 
at a cost which is 2 to 3 times that of ordinaiy tubings, while 
saving 75$ copper r2"f. This is open for further exploration. 
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In the present work, the problem of optimization of the 
fin geometry and inter-fin distance for aluminium has been solved 
for minimum volume and increased heat transfer. Forced flow under 
steady state is assumed. The various constraints on the problem 
are : 

1. Increased heat transfer in comparison to that of a 
bare tube 

2. Flow conditions 

3. 3iot number based on the transverse fin dimension 

4. Boundary layer interference 

5. Number of fins. 

Nusselt relationship has been used for calculating the 
heat transfer coefficient. The radiation effects have not been 
accounted for, as the temperature difference between the fin and 
the environment is small. 

The temperature distribution in the fin has been obtained 
by electrical analogy methods and the heat transfer rate has been 
verified for the best fin. 


1.2 iitziaiubz zzvi::.' 

Carrier and Anderson P 3 "lhave considered the case of 




straight and annular fins of constant thickness and calculated 
the heat transfer, comparing then with bare tubes and thus 
computed efficiencies, 

Baughty j~4~j has obtained a relationship for short fins 
and long fins with and without insulated ends. He has compared 
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the results and plotted curves to compare then with bare tubes. 

Murray and Cano ridge jjTJ have shown that the addition of 
a fin changes the temperature distribution around the tube. This 
change is negligible in slowly moving air. 

Gardner jj3~j obtained equations for heat transferred from 
a finned tube as well as fin efficiency, in the form of Bessel 
functions. He has considered annular fins and spines of constant 
and varying cross-sections in steady state, assuming constant 
surface temperature and one dimensional conduction. The first 
assumption is seriously questionable in the case of heat exchangers 
though it is valid for evaporators and condensers. 

Wilkins £7^ has shown that for the given rate of heat 
transfer at constant base temperature, the optimum fin profile has 
a very sharp point end and a triangular profile uses about 63. 5?= 
less material than the other profiles like straight or parabolic. 

Edward and Chaddock jj8j were the first to account for 
radiation from the finned surface. Radiation can take place between 
tube and fin surfaces, bctxTeen two adjacent fins, fin and enclosure, 
and tube and enclosure. The first two contribute negligible 
amount, and may be neglected, but the last two contribute signi- 
ficant amount to heat transfer. The experiment was for natural 
convection and heat transfer coefficient was assumed to be the 


function of temperature. 


Th. S. Schindt 9j has determined the profile of the fin 
requiring the least material. He considered the contribution of 
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each surface of fin and tube separately to heat transfer, by 
varying the fin size. In order to compare different shapes he 
defines fin effectiveness as the ratio of the actual heat dissipa- 
ted by the fin surface to that which would be dissipated if the 
fin surface wore held at the temperature of the base. He proves 
that the tip of the fin should approach the temperature of the 
surrounding fluid. 

It is of primary interest to recognise the conditions for 
which the finned surface has advantages over the unfinned surface. 
Daughty j~4J compared the heat transfer through the base of three 
fins : (a) short fin, (b) the long fin with insulated end, and 
(c) long fin with heat loss from the end. He concludes that the 
use of too short a fin will violate the assumption of one dimen- 
sional conduction. For heat transfer between air and gases, fins 
are advantageous if the characteristic value (reciprocal of Biot 
modulus) defined by 3q. (l-l) is sufficiently high. 

CV • l/Bi = 2k/(h b), (l-l) 

where k - the thermal conductivity of fin material, 

h - the boat transfer coefficient between fin and 
surrounding fluid, 
b - the thickness of the fin. 

If this value is small (<^<^).02) the use of fins 
results in no increase in heat transfer over that foy a hare tube. 
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Hence the selection, of thinner fins of high conductive netal is 
always preferred. 

Hutcheon and Spalding fio] analysed the problem of using 


an analogue computer for free convection. They assumed that heat 
transfer fron the finned surface was proportional to 5 /4th power 
of difference of temperature . L remarkable feature of their work 
is that they accounted for heat transfer fron the fin tip. Their 
results, compared with a case where heat transfer is proportional^ 
to difference of temperature, have been seen to be more accurate. 

Recently, Dent |llj has solved the case of an annular fin 
for natural convection by the equivalent electrical network method 
using an analogue computer. Dent has accounted for heat transfer 
by radiation also, as suggested by Edward and Chaddock|TVJ by using 
the sane expression for shape factor. He has thus obtained 
fairly good results by developing an iterative scheme for use with 
the computer. 

Danilova and Dyundin |l2] have indicated that heat transfer 
coefficient for finned tube is higher than that for a smooth tube. 
The heat transfer coefficient from a finned tube has been sheva to 
depend on factors such as the surrounding fluid, its temperature, 
heat flux, fin geometry and inter-fin distance. The reduction of 
inter-fin distance to some extent increases the heat transfer 


coefficient. 
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Hutcheon and Spalding M have shewn that there is a 
very little difference in the calculated rate of heat transfer, 
even if the heat transfer coefficient is considered as constant 
instead of variable in free convection, under steady state condi- 
tions. McAdams jl4"*j has said that heat transfer coefficient varies 
from point to point for certain types of fins, hut if average con- 
ductivities and average coefficients are used for calculating the 
temperature gradient and heat transfer rate, the results are in 
excellent agreement with experimental results. 

Eatz et al jl5j put forth an expression for heat transfer 
coefficient for a finned tube as a function of fin outer diameter 
and spacing and later established it experimentally. The expression 
suggested is found to he in excellent agreement with the available 
values for natural flow and condensing conditions. The results 
have been verified for various velocities of air flow. 


1.3 MOTIVATION FOE P33SEFT VICZK 

Optimization of fin geometry and inter-fin distance was 
not considered until 1940. Tin geometry was optimized for maximum 
heat transfer rate maintaining the fin radius constant. Edward and 
Chaddock jjifj optimized the inter-fin distance experimentally using 
annular fins of constant thickness for f ree-convecticn. Their 
results were verified by Tent fllj . 



CHAPTER II 


HEAT TSAHSFERIFD F jQM FINNS? TIDE 


2.1 INTRODUCTION 

The use of finned tubes is advisable only if the 

characteristic Biot number based on transverse fin dimension is 

very small. Biot number should always be less than unity, hence 

the thermal conductivities of fin net d. ' in, 1 3 sho-ild be high. Biot 
number based on transverse fin length : 

h b/k ^<^.0 (2-1 ) 

The material chosen for fin is aluminium which has a 

thermal conductivity of 175 kcal/m-hr-°C. The heat transfer 

2 

coefficient for air is between 10-100 kcal/n -hr- 0 C and for 

2 

water as medium 500-5000 kcal/n -hr- 0 C. Hence the use of aluminium 
fins in the air is justified but for water it is seldom recommended. 
For a fin of thickness 2.5 m the Biot number is, 

Bi = (10) {.0025)/l75 (2-2) 

« 1-0 

As the dominating node of heat transfer fron the tube is 
convection # fins add to the heat transfer rate, when the area of 
heat transfer increases. If the height of fins on a tube is 
comparatively small with respect to tube diameter, the formula 
to calculate the heat transfer rate for a plane wall can also he 



CHAPTER II 


HE‘T T^ANST'SE'ES V C I FINNS"' TU33 


2 1 INTBC3U0TIGM 

The use of finned tubes is advisable only if the 

characteristic Biot number based on transveise fm dimension is 

vexy Gmail Biot number should always be less than unity, hence 

tne taermal co ductivitiea of *ia at i„lo oio_ld be high Biot 
number based on transverse fin length 

h b/k 0 (P-l) 

The material chosen for fin is alumniun w^icn mas a 
thermal conducti-ity of 175 kcal/n-hr-°C The heat transfer 
coefficient for air is between 10-100 kcal/n^-hr-°C ard for 
water as medium 500-5000 kcal/n^-hr-°C Hence the use of alurxmun 
fins in tne air is justified but for water it is sel^o- leoomnended 
For a fin of thickness 2 5 mm tne Biot number is, 

Bi = (10) ( 0025)/175 (2-2) 

« 1 0 

As the dominating mode of heat transfer from the tube is 
convcction^fins aad to tne hoat transfer rate, T"en tie area of 
heat transfer increases If the height of fins on a tune 10 
comparatively small with respect to tube diameter, tne formula 
to calculate t^e heat transfer rate for a plane \all era clso be 
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used for a finned tube. For pin type rectangular fins, heat 
transfer can be calculated by using the fomula for a straight 
rod v/ith appropriate boundary conditions. Using these formulae 
it can be established that the use of fins over unfinned surfaces 
is advantageous if the characteristic value is greater than 5. 

Then, an optimum fin geometry can be predicted [ 13 j • 

Ihen fins are advantageous, the heat transfer can be 
increased by placing then as near one another as practicable. The 
heat transfer coefficient decreases when the boundary layers, 
which occur on the surfaces of the two adjacent fins, mutually 
influence each other. The distance between two fins oust not be 
appreciably smaller than twice the boundary layer thickness. 
Comparison between theory and experiment shows good agreement in 
computed and measured heat transfer coefficients, until the distance 
between fins is less than about l/3rd the 'diameter of tube. But, 
if the distance between the fins is decreased further, the 
difference becomes greater. 

Air flowing along a flat pi abe of 30 cm length idth a 
velocity of 15 m/sec creates a boundary layer, approximately 
2.5 mm thick. A lower velocity ( 1 n/sec) results in larger 

boundary layer, about 12.5 m thick. 


2.2 ASSUMPTIONS IF C/JjCULAPIMt HZ^I 1 TS&1BS32 ??DH FI1E13B TU33 
The following assumptions have been used to calculate 
heat transfer from finned tubes : 
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1. Steady state conditions exist. 

2. Air flowing across fins has uniform temperature . 

3. The temperature of tube surface is uniform. 

4. As the thickness of the fin is snail, as compared to 
fin diameter and inter-fin distance, the temperature 
gradient, normal to fin surface, is negligible. 

5. Ends of fins do not transfer significant amounts of 
energy as compared with other parts of the fin i.e. 
the ends are insulated. 

6. Heat transfer coefficient is the sane all over the 
fin and the tube surfaces. 

7. There is no heat source present in the fin. 

8. Fin material is homogeneous and isotropic and its 
conductivity is constant. 

A few of the assumptions listed above are seriously 
questionable. The hoat transfer coefficient does vary from point 
to point on the fin, hut the use of average conductivity and 
coefficient in the theoretical calculation for heat transfer and 
temperature gradient gives good agreement with experimental 
results [6] . 

The temperature of the tube is measured along its peri- 
phery and then averaged. It justifies assumption 3. 
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The tip area of the fin is negligibly small in comparison 

to the total heat transfer area, and hence it adls less than 2^ 

to the total heat transfer f 8 ] . Thus its effect can safely be 

- J in 

neglected and ends can be taken to be insulated asV ^eeumod fea 5 . 
This will be discussed further while discussing the results. 


2.3 


HFAT T^jjSPgR FRO! I FIN 


Because of relatively high themal conductivity of fin 
metal it has been found possible to consider the temperature uniform 
at any cross section. This means that the conduction is one 

]■ 

Considering Fig. 2-1 , and assuming steady state conditions, 
an energy balance on a differential ring of area 46. yields 



Q s s Q + Q 

^cond. 6 conv. rad. 


(2-3) 


Net conduction into element 


? cond. ' 4*' k ‘ 27rb k I r 


[■ 


d 2 0 do 


dr 


2 + dr. 


dr 


Total convective transfer from the surface of element is 
= 4 JT r h 0 dr 


(2-4) 


(2-5) 


~conv. 

In accounting for the radiative exchanges with the 
differential ring, three separate transfers must be considered ; 

(1) that between the tube and element, 

(2) that between the enclosure and element, 

and (3) that between the element on one fin and on an 


adjacent fin 
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For a small temperature difference between the tube and 
the surroundings, the total radiation energy amounts to a small 
quantity as it is proportional to the difference of fourth power 
of absolute temperature; hence Eq. (2-3) becomes : 


p 

d G 1 
1? ** 


d& 

dr 


2h 

kb 


e . o 


( 2 - 6 ) 


■at the base of the fin the temperature- is the seme as 
that of the tube which leads, along xilth assumption 5, to the 
boundary condition : 


r 


o = e 


v 


and r = r Q 




Equation (2-6) along with the boundary conditions of 


Eq. (2-7) has been solved in terns of modified Bessel functions 



and verified numerically : 

E l( r © f?) M r ^ ) + T 1 C r o > K o ^ ) 

8 = 81 M r = r# 7 h > + h ( r o rip E <- - r i W 1 


(2-8) 


Hence heat flow into the base becomes : 


Q = 2 


where 


_ , . ^ > K 1 (r l% > -%>„n»)h(r 1 l?) 

r * 1 > K o ( r ify > * -V r o ty > 

(2-9) 

^5 - 2h/(k b) (2-10) 

Gardner [_6jj has generalised the above relationship for 


extended surfaces of many configurations. These results are 
useful in the design of extended surfaces and presented in the 
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fora of graph in terns of the effectiveness based on a sinilar fin 
of infinite conductivity. 


2.4 3 i zap ip^-si:?. c:irriciz:r 

The heat transfer coefficient for all surfaces is assumed 
to be the sane. It has been observed that fin spacing affects the 
rate of heat transfer considerably. For any finned tube, the rate 
of heat transfer generally increases with reduced intor-fin dis- 
tance fl7j . However, the volume of fins cannot be increased, as 
there is a limit on inter-fin space; the fins should be spaced in 
such a manner that adjoining fins do not interfere with each other. 

Katz et al |jL?fJ used a modified Nusselt’s relationship to 
predict the heat transfer coefficient for finned tubes with 
Freon-12 condensing around it. The results were found to be in 
agreement with the experimental values, They considered the finned 
tube as horizontal with vertical surfaces attached to it, for natu- 
ral convection. 

Assuming laminar flow xath forced convection, the Husselt’s 
relationship for vertical plate can be used to predict the heat 
transfer coefficient for Reynolds number ranging from 5,090 to 
20,090. 


The Husselt relation for heat transfer coefficient is: 
h = 0.664 k a (Pr) l/3 f (2-11) 
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The characteristic length L can he defined for circular 
fins to be equal to the area of one surface of the fin divided by 
the height of the fin. Thus we get, 


L “ 7T - r^)/|. (r Q - r^. 


( 2 - 12 ) 


For small diameter tubes, as used in the condenser, if 
(r Q - r^) can be approximated as equal to r Q (r^ being small), 

Eq. (2-12) gives the characteristic length as defined by Zatz et al 
The values of heat transfer coefficient obtained theoreti- 
cally using Eq. (2-12) were found to be in agreement with experiments 


for cross-flow 


M 


2.5 FLOT CCrorriONS 

The Reynolds number of flow is defined as : 


Be . UL/t> 


(2-13) 


The characteristic dimension L here is defined as : 

L = (d o + d 1 )/2.0 (2-14) 

For a solid cylinder placed in a stream of air, both the 
outer and inner fin diameters become equal and we obtain the cylinder 
diameter as the characteristic length, 

L = d o 


(2-15) 



CHAPTER HI 


FC2MULATI0H Or T£3 '“FTIMIZ./TION PROBLEM 

3 1 INTBOBbCTION 

Methods of variational calculus can bo used to fi" t o 
optima pointJC^ of a fuiction f(X^) f \i en t' ore are si_plc 
restrictions on the choice of But m practice, rcry reotric- 

tions "ro suporinposcd on tf-c* choice of They linit t_c rn. pe 

of permissible values of ^ , mice the process of finding 
co-’olicate.i Many icthods have been developed to tochle sv.ci 
situations depending pn inrily cn the npturo of the function "rJ 
the root net ms 

The rostnctions inposec. n t r o ch^ce of the iOsi£n 
vector ( A^ , pro celled constraints The constraints g^ , any 

be of the fora of inequality or e^u^lity, sach that, 

<-j (,*,) A o, (3-i) 

0 = 1» n 

The donr in obtained by tho hypv.r surfaces Jefi^ac hv t*_a 
constraints diviaes region into t-o, one of which satisfies ell 
the constraints iodine J by Eq (3-1 \ is cnllod feasible region 
and ovary point m tns regi-n is tre foasiLl^ point The 
optima point is a feasible poirt v ere t’ e function vpluo la 


a ninimu 
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^he direct methods, such as the nethod of feasible 
directions, and indirect nethods are used for non-linear functions. 

Methods of feasible directions consist of finding a 
usable feasible direction and optinun step size for each step of 
the iteration. The techniques of finding a usable feasible 
direction is quite complicated and time consuming. The gradient 
projection nethod is effective for linear constraints. 

Linear programing approach is very useful for solving 
problens where constraints as well as objective function care 
linear. 

The indirect nethods are useful where function and 
constraints are complicated. It includes the U3e of the interior 
or the exterior penalty functions. 

The heat transfer constraint, which uses the Bessel 
functions, makes the problen complicated. The interior penalty 
function is preferred to others, when the initial point can be 
chosen in the feasible region. The subsequent points obtained 
during nininization will then necessarily lie in the feasible 
dcaain, owing to the nature of the penalty function. The uncons- 
traintd ninina, finally converge to the relative minimum. 

3.2 MOTIYjJION FC3 PZNdLTY FUNCTION 

Consider sene constraints on X of function F ( X ). 

/W' 'W\ 


Then, the problen bee ones : 
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Mininize F ( X ) , 
' ^ \,v 1 

subjected to 




(3-2) 


1 ,. 


EU 


The ninine of the function could be ensured by solving 
a set of differential equations : 


<4-^ - - v?<x) 

d cL 'w 


(3-3) 


The function will approach the nininun point, where 

V F - 0 (3-4) 

The effect of constraints is yet to be incorporated in 
Eq. (3-3) and the approach adopted is to add scnething to Eq. (3-3) 
which would bend the trajectory to the feasible region. 

Hence 


doC 


= - VF (^)-Zh. vs ( xj , 


M. = H. ( X ) - 
3 3 w 


g. ( X ) 4C 0, 

3 x jw> ' ^ 7 


r f g. ( x ) 'y- 0. 

^ & j v AAn. / ^ 


(3-5) 
(3-6 ) 
(3-7) 


As 11. is discontinuous, the above set of equations any 

u 

not necessarily give a solution, but this could be a pi at fairs to 
search for a penalty approach and constrained problens could he 
reduced to unconstrained ones as in Eq. (3-5). 
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The penalty function thus constructed is sequentially 
minimized until the optimum is reached, by varying the penalty 
parameter. 

We define the penalty function 

m 

P ( X , r, ) = F ( X ) - r, ^ l.O/g. { X ) (3-8) 

3=1 

where r^ is a penalty parameter. 

The minimization of Eq. (3-8) is sought for a decreasing 
sequence of penalty parameter r^. The penalty function always has 
a minimum in the feasible region, hence the relative minimum is 
forced towards the constrained minimum from the interior region. 

As long as the optimum point, X. is in the feasible 
region, the penalty tem will add a mall quantity to the objective 
function F ( X^ ) , but as some boundary is approached, one of the 
gj ( X^ ) shoots up and the search direction will he changed. If 
r^ is made smaller in each cycle of minimization then, at a point 
X_ the magnitude of penalty tem becomes insignificant compared 
with tho objective function and the point X^ is accepted as the 
minimum of F ( X ). 

The initial value of r^ is chosen, such that, at the 
initial point, both the objective function and the penalty terns 
in Eq. (3-8) have almost the same weight age. Subsequently, r^ 
is decreased in each cycle by some fraction jj20 j . 
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It is observed, that it takes 4-7 cycles to reach 


mmnun. 


ALGGlcITHM 


1. Choose a feasible starting pcint X Q , satisfying the 
constraints 


g . { X ) 0. 

2. Compute initial r^ such that : 


r k = ABS 


F( * 0 )/T. l,0/ gj (Xj 

3 


(3-9) 


(3-10) 




3, Construct the penalty function and obtain the optinun 
point, X, , for P (X, r. ), using sone good technique 

- ** 'ViA- ^ 

iW\i 

for the unconstrained minimization. 

4. If all convergence criteria are satisfied and r^ has 
been changed 3-4 tines, then terminate; 
otherwise, change r^, such that, 


r k«- r k c 
where c 1. 


X 


fVV' 


(3-11) 


(o— 12) 


and repeat fron step 3. 

A flow chart of above scheme appears in Appendix A. 


3.3 UNC0NSTPAIN5D MINIMIZATION 

Many iterative schemes have been developed for the 
unconstrained nininization. The Grid and Sand on methods nay not 
be able to locate exact minima, besides, these methods my need 
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several function evaluations. The zeroth order methods, univariate 
and Powell's ic not require the derivatives of function. Powell's 
method is quadrctically convergent and for the quadratic functions 
it converges in af.ew .steps,, but for the complicated functions it 
requires several cycles of minimization. Methods which use the 
first derivatives cf function, are the first order methods, such 
as, steepest descent, Fletcher Eecves or conjugate gradient 
methods. The steepest descent method is veiy slow because the 
successive moves are perpendicular. Fletcher-Seeves method is 
efficient for the quadratic functions but it takes a considerable 
tine for complicated functions. Newton's procedure and David o;>r 
Fletche r-Povel 1 method are second cider methods. Newton’s method 
needs the second order derivatives cf the function and also inver- 
sion of matrices. The method is found to he strong, but rr>b 
efficient. Sometimes it nay be difficult to evaluate the second 
derivatives of the function. Bavidon-Fletcher-Powell procedure 
does not require the second derivatives of function or nefcrix 
inversion, oven then it is colled second order method, as at 
the optimum point the H-natrix converges to the second cider deri- 
vatives matrix. David on- Fletcher^-Powell method, popularly known 
as D-F-P is quadrat ically convergent. 
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/iliGOEITEU 


1. Assume an initial positive definite symmetric natrixjfs'j , 
sayj^ll ; corresponding to the initial guess vector X , 


, 0 


find the search direction S : 

-vjR 

S = - r H \7P { X ) 

° L OJ v K ,W 


{3-13) 


f\TS/\ 


2. Use a powerful one dimensional minimization nothod to 
compute the optimum step-size 1*. , to obtain the new 


design vector X ^ 

*V\<~v 


(3-14) 


X i + 1 " h ■ 

'•✓W "W 

where qT minimizes F ( X. + (X. S.) one dimensionally. 

*Vv * r. i 

3, If the convergence criterion is satisfied, then terminate; 
otherwise compute a new H-aatrix |JEL + ^ j 

M-W + [\; + H 


as 


(3-15) 


where 


w-orT», s i T j/ ( s i 1 j; - 

-w r v ^ "<W* 

[n. 1 * -f(H. y.)(H. Y. )"*!/( Y. T H. Y. ) 
- V * ( X. +1 ) - V F ( X. ) 

«w\ A^XVV 

4, Compute the new search direction S 


and 


i+1 ' 


(3-16) 

{3-17) 

(3-13) 

(3-1S) 


a^vvv 


and repeat from st.ep 2. 
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The optinun step length, , should be computed accurately 
to ensure that the new H-natrix generated is positive definite, for 
which the method will work efficiently and ensures F^ + ^ ^F.. The 
rounding off errors’ involved in computing 0( accumulate and conta- 
minate the search direction, S^ + ^. Sometimes the process is 
restarted to overcome it |^2J j 

The flow chart is given in Appendix A. 

3„ 4 ONE PESSSIC&JL MINIMIZATION 

Gne dimensional function, F ( o(. ) , can be interpolated by 
a polynomial of any degree and the optimum of the polynomial can. 
be computed. Quadratic interpolation method does not use the 
gradient of function, but it requires function evaluation at three 
points. Fibonacci search and Golden section methods, in general^ 
require many function evaluations and are rarely used in practical 
problems. The direct root method uses the gradient and function 
evaluation at two points on either side of the minimum, and 
interpolate function linearly between these two points. Cubic 
interpolation method can be used if the first derivatives of the 
function are continuous. It uses the gradient and function evalua- 
tion at two points, such that, A o< 3. Each refit mil 

narrow the gap, (B-A) , and ultimately the minimum will be 
reached. This method has been used here and the is computed 
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fron the set : 


,1 


Z = 3 . F { X A ) - F ( Xg ) / (B-A) + 7F ( ) + V F(Zg) 

I * -\*\ 'V\>* t A/y\ 

(3-20) 


r 2 '■ 1/2 

q - r/F(x ) V;-- F (Xg) ( 

I '• '•'/* -V\A J 


(3-21) 


0( “ A + (B-A) ( V F( X A ) + Z+Q)/( V F(X A )+ 7F(X B ) + 2S) 

I vV*' Vv * f Vv\ 

(3-22) 

“d X i+i 3 x i + oC ( S i> ( 3 ~ 23 ) 

'ywv ^VVwr, <VVs -* 

The process of one dinensicnal nininization is repeated 
until convergence is obtained or the nunber of preassigned 
refits is coopleted. 

The flow chart of one dimensional nininization is given 
in Appendix A. 

3.5 COMVEBSBMCB CBITBRIA 

Several convergence criteria can be used to ensure the 
optimum point. The most suitable criteria are those, which are 
applied to the objective function and the design v e ct -e r « These 
ensure that two consecutive minima. do not differ much. These 
convergence criteria can be stated as : 


and 


\ J 

K ) - F (h )] / F < X i + l>! < B 2 

I U- A/W, iVS-w 1 


(3-24) 

(3-25) 
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One of the best criteria used for the termination of 
3-F-p nethod is 


Tv? { \ )1 T I'h. J jVHx.) '/f( x.) 

— **\ 'Nyv L. . 'Vv*. 






-28) 


The gradient at the niniorr, should he zero and it serves 
as one of the- termination criteria for one dimensional cubic 
interpolation. 

\.i"l 1 [v n hJ | £ e 4 < 3 - 27 > 

— rvw\ - -v^ { 

* 

To ensure accuracy of step-size, <<; , one more termination 
criterion the orthogonality test, is used in one-dinensional 
minimization 

$ 

Pr] T F ( X i +l )T /{ ! Sjl ]^( x 1 + 1 )j}!<»» <■*-&>) 

The values of E^, E^ , Eg, E^, Eg have been chosen from 2 
percent to 0..1 percent. 


3-S CCNSTddINTS 

The various constraints enforced, are on heat transfer, 
Reynolds number of flow, Biot number, the number of fins etc. 


3.6-1 HEJiT T&LNSEEE COSSTmUT : 

The rate of heat transfer from finned tube ought to 
be many tines that of a bare tube. For imposing this constraint, 
the heat transfer from the finned tube is written as : 


Q -f>KV 


(3-29) 
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where 3' is the heat transfer from a bare tube under flow condi- 
tions identical to those of the finned tube, 'The value of I 
which is a constant ( l) is usually very high between 4-40, 
hence Eq. (3-29) becomes 

Q ->■ 40 Q’ (3-30) 

if larger value of E is used. 

3.6- 2 PLOT CONDITION : 

The condition of forced convection in laminar flow has 
to be maintained, otherwise the equations used for computing heat 
transfer will become invalid. The characteristic length for calcu- 
lating Reynolds number of flow has been given in Eq. (2-13), To 
have established flow, the velocity of flow can be varied from 
1.5 m/s to 5 n/s and Reynolds number should not exceed 7000. 

Thus : 

B e <C 7C00.0. (3-31) 

3.6- 3 BIOT NBL33R CONSTRAINT : 

The heat transfer area can be increased by increasing the 
tube size or by providing fins on it, While any increase in the 
volume of tube is expensive due to high cost of copper, the use 
of fins of some highly conducting material, such as aluminium is 
advantageous. The limit of the conditions for which fins are 
advantageous is that heat transfer should increase while its 
cost decreases. For this the Biot number, based on fin thickness 
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and outer fin radius, should not exceed unity. Further, the 
accuracy of fin calculations depends on ho w accurately the Biot 
number is computed. 

For putting a limit, the Biot number has been obtained 
for the optimum fin geometry of free convection J 8J and taken 
as a reference, as the Biot number in forced convection will be 
further less : 

h b/k <. .00067 (3-32-) 

h r Q /k < 0.03 (3-33) 

3.6-4 3CnNDJtfff LA2EE IlfTEHFBSBNCE : 

To avoid boundary layer interference on the adjoining 
fins, the inter-fin distance should be slightly smaller than 
twice the boundary layer thickness, but should be appreciably 
large. Comparison between theory and experiment shews that as 
a good agreement between the two, it should be slightly smaller 
than dj/3.0, together with Sq, (3-32), the minimum number of 
fins should be : 

f»-r * *-f r *°' a ***> 

n = 30.0/ \ 0.4 + . 00067(k)/h ) (3-34) 

A 35 J . 

To specify a constraint 

35.5 n 45.5 (3-35) 

The upper limit on the number of fins will check steep fall 
in inter-fin distance. 
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3.6-5 OBVIOUS CSS3I7«;i7S : 

Apart from the so, sene obvious restrictions include 
non-negativity constraints on fin geometry and inter-fin distance : 


d >. d 4 
o l 

(3-36) 

s ^ C 

(3-37) 

b ^ 0 

(3-30) 


The constraints have been normalised. The actual partial 
derivatives for all the constraints and objective function have 
been computed. The finite difference method has been used for 
the heat transfer constraint since Bessel functions are involved 
in its calculations. 3ach function minimum in sequential opti- 
mization is checked for constraints, to make sure of its feasibility. 


3.7 CCiFCTZj. PSCGr^dZ 

The computer programme listing has been given in 
Appendix A. The stepsinvolved are described below : 


1. At each value of design variable, the value of a function, 
f(X) specified below is evaluated : 

W’ 


f Ot ) - 71 (4® - d, > b b/4.0 - r k 2J 1.0/ gj ( S_) 

(3*t3§) 


where g. ( X ) are constraints : 


(3-40) 
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2. The heat transfer coefficient is evaluated before function 
evaluation by using the Nusselt relationship, defined by 

B** (2-11). 

3. The Beynolds number of flow is obtained, by using Eq. (2-13). 

4. The 3essel functions have been confuted fron series, taking 
first 10 dominating terns of series. 

5. For each value of r^, the Davidon-Fletcher-Powell method is 

applied to find the unconstrained minimum. The process is 

repeated till the convergence criteria are satisfied. 

the 

6. The value of is changed till^/relative minimum is obtained, 
which satisfies vector termination criteria. 



CHAPTER IV 


ELECTRICAL ®T.:0?J{ ANALOGY 


4.1 IIITRCSUCTICN 

The analytical solution, discussed in Chapter II for 
combined convective and radiative heat transfer from finned, 
surface is quite difficult as it involves a complex mathematical 
formulation. Often, only the network method can be used, since 
it is quite simple and efficient. It is analysed here, under 
the sane assumptions as listed in Section 2.2. 


4.2 BEAT TRANSFER ~Ry; FIT: 

Consider the fin shown in Fig. 4.1 and a differential 
ring at a radius x-g. Cn writing an energy balance across the 
ring, considering the fin as a cylinder, we obtain : 


W 4 TX * » (®i - ») 3 Ao<5 e ( r / r l>3 - 4 IT k * (® - e 2 I 3 A=>g e (r 2 

( 4 - 1 ) 

Convection from the upper and lower surfaces of the ring : 


Q conv. = 4 7T r B < r 2 " r l>B h c (d “ V® ' 4 ~ 2 > 

As the radial thickness of the node is small, the 
convective heat transfer from both surfaces will be small compa- 
red with conduction into the node and hence the one-dimensional 
conduction heat transfer equation for cylinder can be used. 
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Fig. 4-1 Representation of node on the fin 
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Dent [*ll*j obtained an agreement of 5-15$ with those of experimental 
results using the one-dincnsional heat conduction equation for 
cylinder. 

The radiation from upper and lower surfaces of the ring 

to environment .at temperature & , which is small, can be neglected. 

§ 

Thus, energy balance on the ring : 

foj ~ e ^B (Q - Q 2 ) b _ ( e - 3 g )g 

Iog e (r/r 1 ) B 7(4 TT k b) ~ log e ('r/r 1 ) B /(4 TT k b) “ l.o/ ~4JT\ r B ' ^V^B 


(4-3) 


4.3 :< 5 r:ox; .jljlogy 

The heat transfer equation in the preceding section 
can be equated to a network shown in Fig. 4.2. The application 
of Eirchhoffs law to this network yields : 


v _ v v _ v 

B1 B V 3 V 32 


V V 

B B 

■ ■ " " + tt- 


(4-4) 


^1 ^2 
One can see that Eq. (4-3) and 3q. (4-4) are analogous; 
conparison of both gives : 

1^= log e (r/r JL ) B /(4f]‘k b) (4-5) 

\ 2 = lo g e ( r 2/ r )B^ 4 ^~ ki ^ t 4 - 6 ) 

S II =l/(4 7Tr B (r 2 -r 1 ) B h c ) (4-7) 

Hence it is seen that the fin can be represented by an 


equivalent network as shown in Fig, 4,3 




Fig. 4.2 : Network r-e-pre: 



S'ig. 4.3 : Network retire sei 



ion of the fin 



34 


4.4 rrrrjra: sche:s res solution cr ozricrff 

To establish the temperature distribution in the fin, 
for an assuned value of the temperature distribution, heat 
transfer is computed by calculating the network resistances of 
the circuit. Cnee the total heat transfer (or current flow in 
circuit) is computed, actual temperature distribution through the 
network can be obtained. If it does not agree with the assuned 
value, the process is restarted with computed temperature distri- 
bution as the initial guess. The cycle is repeated until 
agreement is reached between the assuned and the computed tempe- 
rature distribution. The convergence to the correct solution 
is veiy rapid. 

Uhen convergence is satisfied, the problen for total 
heat transferred by the fin has been solved and hence fin effi- 
ciency could be computed. 

The number of nodes depends on the problem. As in other 
numerical techniques, the more the nodes used, the nigher the 
accuracy but the more the labour involved in solution of the 
problem. Usually 5 to 8 nodes can be used as a compromise 
between accuracy and computer tine involved. 

A computer programme written for above iterative 
scheme appears in Appendix B. 



CHAPTEIt V 


SESULTS MID DISCUSS 101® 


5.1 INTRODUCTION 

The optimum fin geometry and inter-fin distance have been 
obtained using optimization techniques for different temperature 
differences between the tube and surroundings. The results are 
presented in Table 1. Since the available wind tunnel for experi- 
mentation does not permit velocity of £107; more than 4 m/s at the 

test section, and for laminar flow Reynolds number is to be maintained 

> 

7000, all calculations have been carried out for a velocity of flow 
of 1.7 m/s. 

The electrical network analogy has been used to obtain 
the radial temperature drop in the fin and verify the heat transfer 
rate. 

5.2 :-IZ jg TDu^CSrZTi CCEFFICLIaT 

The heat transfer coefficient was computed using Dusselt 
relationship for a vertical plate. The characteristic length has 
been defined in Eq. (2-12). Since the heat transfer rate from the 
tube is small (40 to 50 tines!) in comparison to that of fins, the 
heat transfer coefficient for the tube was also assumed the same, , 
Figure 5-1 shows the variation of Musselt number with 
l/(Se Pr) for Prandtl number of 0.72 as calculated by using 



NUS SELT No. 



VO/mePr) 

Fiql ;:5-1 Vo r iati^r* of Nussett number with l-D/fRePrifor 4hef#l 



37 


tho Eq (2~ll) The results aro found to bo in agreement with 
those of constant rail temperature vortical plate ^19 ( 

Table 2 shows tho effect of the inter-fin distance on 
tho heat transfer coefficient The heat transfer c^efficieit 
irereascs to some extent, if the inter- fin distance is decrease,^} 
t rise of about 6$ m the heat transfer coefficient is observed for 
a decrease of 67$ m the inter-fin dist jaco If the inter-fin 
distance is loss than d^/l2 0 then, on increase "in heat transfer 
coefficient is observed, oien if mtoi-fm distanct is incre^sec. 

This establishes the effect of boundary layor interference on iiite-- 
fin distance and the heat transfer coefficient J" 2^ 

6 3 TSijpEP^JCnB DISTRIBUTION IN THG FIN 

T utcheon and Spalding | 10j assn td the vonotion o* 11 o 

temoe rature as 5/4th power for tao network method The radial 

temperature distribution obtained by tho nuroricol method and ti.o 

✓ 

network Liothods are jresontod in '’’aiLe 3 and plotted in Fig S- 9 

An excellent agreement, v/ithm 2 peicent, 1b foun bet roen 
tho two methods 


b 4 OPTIMUM SOLUTION 

Tables 1 and 2 contain the various values computod by 
the optimization programme It is observed th^t the process 
converges to the minimum m 3 to 6 cycles and takes about 8 to 
9 minutes, tho compilation time of the programme being 4 minutes 



TEMPERATURE OF NODE TO TUBE TEMPERATURE 



-4 ; J i i _l 

2 3 4 5s 

(r/r t ) 

RATIO OF NODE RADIUS TO TUBE RADIUS 


-L 

8 


Fig, 5~2 Comparison of the radial temperature drop 

in the fin by the network and the actual methoi 
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The execution tine depends on the initial guess, -which was the sane 
for all temperature differences between the fins and surroundings. 

J-he penalty function value has been reduced by 80—85 percent and 
tne volume of the fins decreases by 40—50$. The heat transfer rate 
in comparison to that from bare tube under sane condition of flow is 
40-50 tines no re. 

5 . 5 HZ-'J TR'JS'SraHr.ZD FPCII TKZ Til CT T. - ~ SIT 

The tip of the fin was assumed to be insulated for 
deriving boundary conditions for Eq. (2-6). Now accounting for the 
tip area of the fin; for the fin of maximum thickness : 

Ratio of 'tip area /fin area = Jf (l0.83)(0. 127)/ TT (l3.88 2 -l.£2 2 )/2.0 + 

77(10.88) (0.127) , 

= 1.93$ (5-1) 

which can bo neglected in comparison to the total surface area. 

5.6 H3ZT TR'JSBEZEZZD FROM THE FINNED TUBS 

Figure. 5-3 shows the variation of the boat transfer with the 
temperature difference. The network method has been used to compute 
the heat transfer rate for the best fin. The deviation of 0.13$ to 
1.92$ is observed in the heat transfer rate computed by two nethedb. 

The agreement is excellent, even when the equation for one dimensional 
heat transfer from a cylinder was used. This justifies use of the 
network method in comparison to the actual method which involves too 
much of computer time, while the network method is veiy efficient 


and fast 



RATIO OF HEAT TRANSFERRED BY FINNED TUBE TO 
BARE TUBE AT REF TEMP 10 3 C 


REFERENCE TEMP 10°C 
• NETWORK METHOD 


METHOD 


exA^T 


( 0,/e ) 

RATIO OF THE TEMPERATURE C?F FINNED TUBE TO THE REE TEMP. ?0*C 


Pig. 5-3 Qqmpar isch of the heat transfer rate by the network 
andwlKe actual methods 
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?A3I3 


OPTIMUM PIN 33r.‘3T3Sr AND II'TXSR-EIIT DISTANCE 



T 
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$ 

IMIjEZ 
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IITEX-JUT £ 
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ul j «oj 
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43 

13.37 

0,83 

r* oo 

v> • Jo 

531.0 

o 

6750 

43 

10.27 

0.503 

0.127 

402.2 

200.5 

6073 

44 
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0.53 

0.127 

525.2 
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TA3LS 2 

CPTI'ilZATIGN RESULTS 


ttjycla 
Tem- Jpf un- 
pera-|: ons- 
ture grained 
0 foinimi- 


Penalty | Penalty f 7olun€$Tc . of| Heat |lnt s r— | 3s at 
jaraneter $ function | 3 | f in^t trams— jjfin § trans 

| | Cn I | ferred f dis- 5 fer 


trans- ffin § trans- *efie- ^ muter 
ferred * dis- $ fer $3tr>re-£ t~.n§ 
kcal/ % -i| tance * coeff . ness * rts^ 

| ca | l£cal/jhr* } 


10 

1 



1.89 

3 

-4 

ID 

7.1 

3< 

IO" 2 


i t 

39 

it 

225.3 

•j , 573 

L. “ C S 

1.285 

10“ 3 

t 

— - ■■■■-" —■ " 


2 

9.5 

10 -5 

4.6 

io” 2 


40 

255.8 

3.57^ 

1.281 

t "~r3 

i. -J 




3 

4.74 

io" 5 

2.6 

io -2 


43 

243.8 

0.560 

1.30 

10 r3 




4 

2.0 

10" 5 

2.2 

io -2 


44 

220.0 

0.563 

1.33 

io- 3 




5 

8.5 

!0" 6 

1.39 

io -2 

391.8 

43 

23-3.2 

0.830 

1.345 

io- 3 

70.4 

9 

0.3 

1 

1.08 

_ 4 

10 



742.0 

41 


0.568 

1.356 

13~ 3 




2 

9.42 

10" 5 



668.0 

42 

225.5 

0.566 

1.37 

10~ 3 




3 

4.71 

10" 5 



570.0 

43 

222.5 

0.565 

1.38 

io- 3 




4 

3. OS 

10" 5 



416.0 

44 

213.5 

0.567 

1.385 

jlO 




5 

1.99 

10" 5 



412.0 

44 

211.0 

0.568 

1.388 

10~^ 




6 

1.29 

10" 5 

1.85 

io" 2 

402.2 

43 

200.5 

0.576 

1.33 

10 -3 

70.2 

9 

7.2 

1 

1.9 

io" 4 

5.7 

io- 2 


40 

212.0 

0.60 

1.285 

1c -8 




2 

9.3 

io- 5 

4.0 

w- 2 


33 

211.8 

0.573 

1.233 

o 

10“=” 




3 

4.7 

10" 5 

3.56 

io- 2 


44 

220.4 

0.564 

1.32 

« r~Z 

JlSj > 




4 

3.04 

10‘ 5 

2.98 

io -2 

525.2 

44 

221.0 

0.56 

1.35 

io- 3 

78.0 
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TABLE 3 


RADIAL TEJPILArjLB DEO? 


, — -ITT-Til 




EEFOiKI 


TIDE AID SBISOUIOING 


RADIUS C? TU3E 


0,61 cn 


RADIUS 0? FIN 


5.40 m 



Actual Metho; 


Network Metis 


3.688 


0.658 


3.656 


0.639 


0.631 


0.632 




Tij3L3 4 


CCiirAilBCN CF F2S 3342 T2&I13732 2413 PCS 
OPTIMUM FILT 33 S SISC &L 1 T5L5G3£jlX3 3I7F3.“3i 1 03 C' 


Te^ip./Eef. Tenp. j 
6/9 i | 

r Actual Method 

i q/qj 

i 

nun 


0.3?. 

1.041 

1.93S 

1.131 

0.94 

1.064 

1.361 

1.131 

0.30 

1.088 

1.083 

“HOC 
-L . J.O fw 

0.98 

1.110 

1.105 

1.132 

1.00 

1.132 

1.135 

1.132 

1.02 

1.155 

1.151 

1.132 

1.04 

1.179 

1.172 

1.132 

1.06 

1.20 

1.198 

1.132 

1.08 

1.223 

1.22 

1.132 

' 1.10 

1.245 

1.242 

1.132 




CHAP TER. VI 


CCNCLU5ICWS AID SCCJ3 73 Y3Z322 I7G3Z 

3.1 CCITCLZBICS 

Ac attempt a as been made to apply tie optimisation techniques 
to tf_.e problems ox beat transfer from fins* In tbs absence of the 
results for forced collection* no comparison with available information 
was possible, nowever * looking to tbe present practice of the use of 
thin fins placed closed to each other in evaporator and heat exchanger, 
our results seem to be quite satisfactory, 3ven then an experiment o.l 
verification is very important. Owing to the nature of the penalty 
function, the minimum obtained cannot be accepted as a relative minimum 
but it always lies very close to relative minima, 

The minimum obtained cannot be taken as an absolute minimum, 
to locate the absolute minimum, the initial guess is to be altered, to 
obtain several relative minima, if they exist. 

Verifications by the network analogy method for the best fin 
has been found to be in agreement of between 9.2^ to 2 JS, which is 
quite satisfactory. The network analogy method can be applied to 
compute the heat transfer rate and the temperature drop in the fin as 
the actual method requires much computer time. 

The computer programmes are quite general and can be used for 
other problems with a little of modification. 

3.2 scope ?C5. Fangs tg 

The experimental verification on the available wind tunnel 
which can produce a flow of 1,5 m/s to 5 n/s can be taken as an 



extension to this p rob lea. Tvro 4C cn long tubes provided vith t*.e 
aluminium fins cen be used. The radial temperature of the fin ana 
tube can be measured by the thermocouples. The velocity of fio*? 
measured by an orificeaeter is varied by an air blouer. 3iectnc or 
steam heating is available for experimentation. 

The effects of radiation have been neglected as the tempe- 
rature difference ?ras small. Optimization techniques may be usen ior 
high temperature differences by, accounting for the radiation. 
shape factor can be obtained by solving the differential equation fbl 
with the help of some numerical technique. 
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